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ABSTRACT 



X 

We discuss the possibility of testing the principle of equivalence with solar 
neutrinos. If there exists a violation of the equivalence principle quarks and 
leptons with different flavors may not universally couple with gravity. The 
method we discuss employs a quantum mechanical phenomenon of neutrino 
oscillation to probe into the non-universality of the gravitational couplings 
of neutrinos. We develop an appropriate formalism to deal with neutrino 
propagation under the weak gravitational fields of the sun in the presence 
of the flavor mixing. We point out that solar neutrino observation by the 
next generation water Cherenkov detectors can improve the existing bound 
on violation of the equivalence principle by 3-4 orders of magnitude if the 
nonadiabatic Mikheyev-Smirnov-Wolfenstein mechanism is the solution to 
the solar neutrino problem. 
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Experimental test of the principle of equivalence, one of the fundamental 
building block of Einstein's general theory of relativity, now has history of 
more than a century. It was Eotvos JTJ who carried out the celebrated torsion 
balance experiment that elevated the experimental test of the equivalence 
principle to the realm of precise measurements. Eotvos and his collaborators 
obtained the bound t] < 3 x 10~ 9 for the violation of the equivalence of inertial 
and gravitational masses. Forty years later the Eotvos-type torsion balance 
experiment was substantially improved by Dicke and his collaborators [[J. 
They noticed that sun's gravitational field affects this type of experiment 
by producing torque of 24 hours period in the presence of violation of the 
equivalence principle. With their extensive efforts in removing systematic 
uncertainties they achieved the accuracy r]<3 x 10~ u , an impressive im- 
provement over two orders of magnitude over Eotvos'. Using the similar 
apparatus with Dicke et al.'s Braginski and Panov |3[] reported the bound 
77 < 0.9 x 10~ 12 , the most stringent one to date. 

In this paper we point out that by using the sun as the source of grav- 
ity and at the same time as the source of neutrino beam one can obtain 
much more stringent constraint for violation of the equivalence principle. 
We show that the optimal sensitivity one would achieve by the solar neu- 
trino observation by the next-generation water Cherenkov detectors is | A/| ~ 
10 -15 — 10 -16 , an improvement over 3-4 orders of magnitude than the Eotvos- 
Dicke-type experiment. Here Af is a measure for violation of the equivalence 
principle for neutrinos to be defined later. 

A mild bound |A/|<10 -3 for violation of the equivalence principle for 
neutrinos have been derived [|J by using the small difference between arrival 
times of photons and neutrinos from SN1987A. The best bound deduced 
so far for microscopic objects is from the neutron free fall refractometry 
experiments which led to \rj\ < 3 x 10~ 4 ||. 

Our discussion in this paper heavily relies on the basic observation by 
Gasperini M . If Einstein's equivalence principle is violated gravity may not 
universally couple with neutrinos with different flavors. He pointed out that 
if this occurs neutrino oscillation similar to that of flavor mixing takes place 
and the effect can be detectable by experiments. After his observation sev- 
eral attempts have been made to sharpen his proposal and to examine the 
attainable sensitivity for violation of the equivalence principle. 

Halprin and Leung [[?[] observed that one can use solar neutrinos to per- 
form a sensitive test of the equivalence principle and they gave an order- 
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of-magnitude estimation of the sensitivity, |A/| ~ 10~ 14 . Iida, Minakata 
and Yasuda || discussed the possibility of doing the similar test by the 
long-baseline accelerator neutrino experiments and obtained |A/| ~ 10~ 14 as 
expected sensitivity. It is accidental that two methods yield the same order 
of magnitude for the sensitivity since the ingredients involved in these esti- 
mations are entirely different with each other. Some related considerations 
were given in Ref. @. 

If one want to use solar neutrinos as an experimental mean to test the 
equivalence principle one has to decide his attitude to the solar neutrino 
problem | ID| . One might say that it is due to the fact that our ability of 



modeling the sun is so poor and the problem may be solved by inventing 
an appropriate non-standard solar model. Others may argue that it implies 
a new physics of neutrinos beyond the standard electroweak gauge theory. 
Someone would be more brave and may propose that it could be explained 
by the effect of violation of the equivalence principle itself, as has been done 
in Ref. @. 

In this paper we make an assumption that the Mikheyev-Smirnov-Wolfenstein 
(MSW) mechanism |Tl| is the cause of the solar neutrino deficit. It makes 
our analysis model-dependent but it is inevitable if the solution to the solar 
neutrino problem involves the physics beyond the standard model of elec- 
troweak interactions. In view of the progress in solar neutrino observation 
done and to be done by the present and the future detectors it would not 
be too optimistic to speculate that we would finally find the answer to the 
solar neutrino puzzle. The point we want to stress is that once we have 
the solution one can perform similar analyses as ours to test the equivalence 
principle and the present work may serve as a prototype for them. 

In this paper we will make a substantial improvement in the treatment 
of possible experimental test of the equivalence principle by solar neutrino 
observation. A part of this work has been presented in Ref. [0 though in a 
very preliminary stage. 

The points we are going to make in this paper are: 

(1) We shall develop a new formalism which describes neutrino propagation 
in an arbitrary spherically symmetric metric of a perfect-fluid star. It enables 
us to calculate neutrino flavor conversion probabilities that are accurate to 
order G, the Newton constant. 

(2) Our analysis of the effects of violation of the equivalence principle relies 
solely on the spectral shape of the 8 B neutrinos. We fully treat the effects of 
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the MSW transformation due to the neutrino flavor mixing as well as that of 
the non-universal gravitational couplings. We discuss how one can separate 
the latter effect from the former's and we estimate the sensitivity by fully 
taking into account of this mixture of two different flavor transformations. 

First let us recall the definition of the equivalence principle. According to 
the text book by Misner, Thorn and Wheeler |D|], the equivalence principle 
is stated as "In any and every local Lorentz frame, anywhere and anytime in 
the universe, all the (non-gravitational) laws of physics must take on their 
familiar special-relativistic forms." It follows from the principle that if the 
gravitational couplings of different flavor neutrinos are different the equiva- 
lence principle is violated. This is because in a local Lorentz frame where 
electron neutrinos obey the familiar special-relativistic wave equation muon 
neutrinos do not. 

In this paper we work with two neutrino flavors and denote them as 
electron and muon neutrinos. We consider the Lagrangian 

i=l,2 

- e { G i) m iVMiV M i 
i=l,2 

+ {interactions with electroweak gauge fields) , ( 1 ) 

where e afl (Gi)(i =1,2) are the vierbein fields of some background metric and 
e(Gi) = dete^Gj). The Newton constants Gi(i = 1,2) are allowed to be 
different for neutrino basis vqi and i*G2 ■ The last term in (JiJ) is assumed to be 
written by flavor (or gauge) eigenstate. We take the most general ansatz that 
the flavor, the gravity and the mass eigenstates are all different with each 
other. The latter two are related with the flavor eigenstate in the following 
way: 



cos 6q sin 9 G 

— sin Q G cos Qq 
cos 9 m sm9 M 

— sin 9 m cos 9 m . 



VM2 



(2) 



The mixing angle 9m in @ is the usual flavor mixing angle which plays an 
important role in the MSW mechanism. 
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As a measure for violation of the equivalence principle we define 



We believe that Af defined in (|3]) does have a right correspondence with r\ 
which is defined in a similar fashion as @ by replacing Gi by (M/m)i as a 
measure for the in-equivalence of inertial (m) to gravitational (M) masses. 
The equivalence of i] and Af comes closer to be true if we are allowed to 
boldly interpret the Eotvos-Dicke-type experiments as measuring the non- 
universality of gravitational couplings of u- and d-quarks. 

We shall derive the neutrino evolution equation under the influence of a 
weak gravitational field of the sun. For generality we work with the following 
most general static and spherically symmetric metric, 

ds 2 = g tt (r)dt 2 + g rr (r)dr 2 - r 2 (d6 2 + sin 2 6d(f) 2 ), (4) 

taking the center of the sun as the origin of the coordinate. We make the 
ansatz that the energy-momentum tensor of the solar matter is given by the 
form of perfect fluid, 

T lu , = Pg lu , + (P + p)U /t U v , (5) 

which should give an excellent approximation to the sun. In (|5|) p and P 
denotes the energy density and the pressure, respectively, and is the 
usual velocity four vector. The metric g rr (r) and <?#(r) can be obtained by 
solving the Einstein equation and are given up to the order of G by |T3| 



g tt {r) = 1 - 2G £ 



M(r') 



dr', (6) 



where 



Mir) 

-g„(r) = 1 + 2G=-H (7) 



M(r)= / dr'A7rr ,2 p(r') (8) 



o 



and it can be interpreted as the total mass contained in the volume with 
radius r for Newtonian stars. In case of the sun, we can neglect the pressure 
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term in equation (|6]) which is very small compared to the term A4(r)/r 2 . In 
fact one can show that the pressure term is of order G 2 by noticing the equa- 
tion of hydrostatic equilibrium, r 2 [dP(r)/dr] = —GAi(r)p. This equation 
holds in most of the stars for which the Newtonian approximation is valid. 
We cannot neglect the pressure term in such stars as white dwarfs or neutron 
stars. We will neglect the pressure term in the following discussions. 

Here we give some clarifying comments on the theoretical basis and the 
assumptions we make in doing pilot computations for testing the equivalence 
principle by using solar neutrinos. First of all we ignore the effects of gravita- 
tional field of the earth, where the neutrino detector is located, by assuming 
that the sun's gravitational field gives a dominant effect. Strictly speaking we 
assume in our treatment that the detector is located at far apart from the sun 
where the metric is approximately Minkowskian. We also ignore the possible 
MSW transformation in the earth which may not give a good approximation 
with the large-mixing-angle solution of the solar neutrino problem. 

One of the most important issue which remains to be explored is the 
problem of possible frame- dependence in the result. We shall make some 
comments on it while it will be addressed in more detail in Ref. fHj| . Since 
the general coordinate invariance is broken in the theory defined by (1) there 
is no guarantee that the bound we obtain is independent of the frame we 
choose. Namely, we will do our computation by taking the inertial frame of 
the sun, but the result of neutrino flavor transformation might be different 
if we do the computation in the inertial frame of the earth even if we keep 
ignoring the effect of earth's gravitational field. 

We argue, however, that the corrections due to the frame-dependence is 
of the order of Af, and is negligibly small if we restrict ourselves into, for 
example, |A/| < 10~ 10 . To this goal let us assume that by taking limit of 
Af — > the theory defined by (1) approaches to general relativity coupled 
with matter-gauge fields. We are to deal with weak gravitational fields of the 
sun and confine ourselves into the perturbation expansion up to first order 
in G. In this weak-field approximation the above assumption should be true. 
Then in the limit of Af — > there should be no frame dependence as it is the 
general relativity. In fact, it is a Newtonian gravity with special relativity. 
We conclude that within the weak-field approximation the frame dependence 
comes in as the first-order in Af and would leave only negligible effects in 
observables. In particular the coefficient of Af is frame independent. 

Let us consider the case that the electron neutrinos are produced at the 
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origin and propagate toward the solar surface. In this case we can choose the 
coordinate with 9 = <fi = to describe the trajectory of neutrino propagation. 
Under this choice of the coordinate, the Dirac equation is given by 



i(e u t )- L 7o d t + i(e')-Vdr-m 



0. 



(9) 



where we note that the vierbeins e\ and are given by the metric as 



(el) 2 = g u and (el) 2 = g rr . Here we have ignored the derivative terms 
in spin connections because it is of the order of ~ 1/ER Q ~ 10~ 22 for E=10 
MeV and is negligibly small. 

By taking into account only the positive energy component of the Dirac 
equation and choosing outgoing wave, we get the evolution equation in the 
weak gravitational field of the sun for one generation of neutrino, 



where 



d 

i—v = E 
dr 



-G 



m 
2E 2 



~M(r) 
— — - 


r°° Mir) ' 




r 





(10) 



(11) 



It is easy to generalize this equation to the two flavor case described 
by the Lagrangian ([[]). The evolution equation in the sun for two neutrino 
flavors is given by 





= H 




- *v 




- *v 



where 

H = 
with 



.d_ 

dr 



V2G F N e + A M sin 2 9 M + A G sin : 
|(A M sin2e M + A G sin 26 G ) 



(12) 



G 



|(A M sin20 
A 



. ^ _.';,/ + Ag sin IQq 
■m cos 2 9m + A G cos 2 9, 



A M =Am 2 /2E, 
A G =Af<j>{r)E. 



(13) 

(14) 
(15) 



where Am 2 = m 2 — m\. Note the difference in energy dependences of Am 
and Aq. This makes it possible to distinguish the effect of the violation of 
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equivalence principle from that of the MSW effect by examining modulation 
of the neutrino energy spectrum. 

With one more ingredient having added to the MSW mechanism the 
resultant equation of neutrino transformation fll2|) is far richer than the usual 
MSWs. We only make a few remarks on some distinctive features of it. The 
resonance condition is obtained by equating the two diagonal terms in the 
Hamiltonian matrix (|T3|) , 

V2G F N e - A M cos 26 m - A G cos2fl G = (16) 

From ([TBI) one observes several notable features: (1) The resonance condition 
can be satisfied even for Am 2 = if Af < 0, the gravitational MSW effect. 

(2) Resonance can occur at outside the sun where N e = if Am 2 ■ Af > 0. 
This should not come as a surprise because <j>(r) plays the role analogous to 
the position-dependent electron density in the MSW effect. 

(3) The mathematical structure of the equation (|T^) is quite different from 
that of the flavor MSW mechanism by having r-dependent term in the off- 
diagonal element in the Hamiltonian flT3]). This makes the usual Landau- 
Zener analysis difficult, or at least highly nontrivial. 

Notice that not only the condition (|T^) but also the adiabaticity condi- 
tion is needed for sufficient resonant conversion of neutrinos. We can obtain 
the adiabaticity condition by demanding that the oscillation length at a res- 
onance point is much shorter than the width of the resonance region. For 
resonance in the interior of the sun the adiabaticity condition reads, 



h (AAfsin20Af + A G sin20 G ) 



» 1, (17) 



2vr (A m cos 26 m + A G cos 26 G ) 
where 

A G = A G + h(dA G /dr). (18) 

To express the adiabaticity condition in the form of (|TTD we have approxi- 
mated the electron number density N e in the sun as 

N e (r) = N exp(-r/h) (r > O.lr ). (19) 

with h ~ t©/10. 

Now we will discuss how sensitively violation of the equivalence principle 
affect the solar neutrino observation. To compute P{v e — > v x \ E u ), the prob- 
ability of observing v x (x = e and r) at the earth, we employ the following 
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procedure. We numerically integrate the evolution equation ([12]) from r = 
to r = lOr where the gravity effect is negligible, |A(j| << Am, for regions 
of parameters of interest to be specified below. We ignore the term A G at 
r > lOr and evaluate P{y e — > v x \ E u ) by 

P{y e -> u e ;E u ) = (l - ^sin 2 2^ M ) |AJ 2 + ~ sin 2 26 M |AJ 2 

- sin20 M cos20 M Re [A„ 6 A*J (20) 
P(i/ e = 1 - P(i/ e -> i/ e ; £„) (21) 

where and A„ M denote the probability amplitudes of neutrinos being v e 
and at r = 10r Q , respectively. 

We calculate the electron energy spectrum f(E e ) to be measured at the 
Super-Kamiokande by the following equations 



f(E e ) = £ [dE u e(E e )[da„ xe (E u ,E e )/dE e 



x P(V e -> ^; ^)[d0 8B (^)/^], (22) 

where da Uxe (E u , E e ) / dE e (x = e,fi) are the neutrino-electron elastic scatter- 
ing cross section and d(j)8^(E u ) / dE u is the differential flux of 8 B neutrino 
calculated by the standard solar model (SSM) [|15|] . 

For the Sudbury Neutrino Observatory (SNO) detector we use a similar 
expression as (|22"D , the one in which one restricts the summation over flavors 
into electron and replaces the neutrino-electron scattering cross section by 
the electron energy distribution for a given neutrino energy, dD(E u , E e ) / dE e , 
for the reaction v e + d^p + p + e~ calculated by Nozawa |TB| . It includes 
the effects of (a) the Fermi motions of nucleons in the target deuterons, 
(b) the final state three-body kinematics, and (c) the final state Coulomb 
interactions. We have used very simplified form of the detection efficiency 
e(E e ) as 

e{E e )=8{E e -E th ), (23) 

with Eth = 5 MeV for both of the detectors. Alternative forms of (p3|) 
may not alter our conclusion because the effect of gravity exists in high 
energy part of the neutrino spectrum. We have taken into account the ex- 



pected energy resolution 14%/y P/lOMeV for Super-Kamiokande |L7] and 
10%/JP/lOMeV for SNO g 
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As a prototype analysis in this paper we pick up the two sets of parameters 
of flavor mixing as representatives of (i) the nonadiabatic and (ii) the large- 
mixing-angle MSW solutions: (i) Am 2 = 6 x lCT 6 eV 2 , sm 2 26 M = 0.01 and 
(ii) Am 2 = 2 x 10" 5 eV 2 , sin 2 26 M = 0.7. 

We do not make any attempt to average over the production points of 
8 B neutrinos but assume that they are produced at the center of the sun. 
Given the parameter of interest the resonance points are well outside the 
production region of 8 B neutrino and therefore our treatment should provide 
a reasonably good approximation. 

To have a feeling on how sensitively violation of the equivalence principle 
affect the MSW-modulated neutrino spectrum we perform a pilot computa- 
tion by taking Qq and 9m equal for the case of nonadiabatic MSW param- 
eters (i). In Figs. 1 and 2 we show (a) recoil electron energy spectra, and 
(b) electron energy spectrum divided by the one predicted by SSM. Fig. 1 
is for Super-Kamiokande and Fig. 2 is for SNO. We have done the compu- 
tation with five parameters, Af = 0, ±5.0 x 10~ 16 , and ±1.0 x 10 -15 . In 
Figs. 1(b) and 2(b) the ratios are normalized at 9 MeV of electron energy. 
Super-Kamiokande is expected fl7j to have accuracy of 2-3% for the ratio of 
modulated spectrum to that of SSM after its 3 years operation depending 
on the energy region and the MSW parameters . Therefore, from Fig. 1 we 
expect that the effect of |A/| of the order of 10~ 15 can be distinguished from 
the pure nonadiabatic MSW solution. 

The similar computation can be done with the parameters (ii) of the 
large-mixing-angle MSW solution. But we do not show the result here for 
reasons which will be explained after the full analysis with the large-mixing 
angle MSW solution. 

A technical comment is in order. Because of the limitation of the com- 
puter time we cut off the lower end of integration over the neutrino energy 
in (^) at E v — 5 MeV. Since we take into account the energy resolution of 
about 1 MeV (Super-Kamiokande) and 0.7 MeV (SNO) at electron energy 
of 5 MeV we would have to include the contribution from neutrino energy 
E u >3 MeV. It thus provides an unconventional way of smoothing out the 
effect of the sharp (step function) cut-off ( p3| ) of the detection efficiency. We 
believe that it does not gives rise to any serious effects to our analysis because 
we use the ratio of the energy spectrum to that of SSM. 

Now we enter into a full analysis with the nonadiabatic MSW solution of 
the solar neutrino problem. We aim at identifying the region of parameters 
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Af and sin 26g in which the spectral shape of 8 B neutrinos can be distin- 
guished with 90% confidence level from that predicted by the pure MSW 
solution with the particular set of parameters (i) and (ii). At the present 
stage it is of course impossible to pinpoint the MSW parameters even if we 
assume that it is the solution to the solar neutrino problem. The present work 
with the particular set of parameters, therefore, is only meant to present a 
proto-typical analysis toward the more complete ones. 

We define, as a quantitative measure for the deviation of the electron 
energy spectrum, 



X 



mm 



N 

E 

i=l 



a; i (A/ = 0)-aa; i (A/,sin2g G ) 
o"i(st) 



+ 



' 1-a V 



(24) 



following the procedure of Ref. [ITS]. We divide the energy range into N 
= 19 bins and calculate the count rate by integrating the electron energy 
spectrum in each bin. Xi in ( p4|) implies the ratio of the count rate calculated 
in this way to the one calculated in the same way but using SSM without 
any effects of flavor-gravitational transformation, erj(st) and <Tj(sys) denote 
the statistical and the systematic errors, respectively. The symbol min in 



front of the right-hand-side of fl2"4 ) is to take minimization by varying the 
parameter a. The measure \ 2 is thereby sensitive only to the shape of the 
energy spectra within the uncertainty of absolute normalization represented 
by the systematic error. 

In our computation we treat the statistical error in the following way. We 
suppose that the Super-Kamiokande (SNO) detector observed 20000 (6000) 
and 40000 (12000) solar neutrino events, which roughly corresponds to its 
operation over 2 and 4 years, respectively. Then, it is rather straightforward 
to evaluate the statistical error by computing number of events in each bin. 
We combine the errors of the numerator and of the denominator by assuming 
that they are independent Gaussian errors. Lacking detailed knowledge of 
the detector performance we simply assume, as systematic errors of both the 
detectors, <7j(sys) = cr(sys) = 6%, the systematic error of the Kamiokande-II 
experiment |T9|| . 

We emphasize that our treatment of errors does not have any particular 
significance and should be taken only as tentative. In particular we did not 
attempt to make a detailed comparison between the expected sensitivities 
between the two detectors. 
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In Fig. 3 we present the results of the analysis for Super-Kamiokande in 
cases of (a) Af > and (b) Af < 0. Plotted in Fig. 3 are the regions of 
parameters on Af — sin 26q plane which are sensitive at 90% confidence level 
to the effect of violation of equivalence principle. Notice that the relative sign 
of $g with respect to 6 m is physically meaningful. The open circles and the 
asterisks are for cases with 40000 and 20000 events, respectively. We observe 
that the sensitive region goes down to | A/| = 10~ 15 — 10~ 16 depending upon 
sin. 29c- It is notable that a factor of 2 better statistics achieved by taking 
the data for 2 more years does not significantly improve the sensitivity. It is 
largely due to a sharp growth of chi-square near the edge of the sensitivity 
region. 

The same analysis is repeated for the SNO detection and the results are 
presented in Fig. 4. The notations are exactly the same as in Fig. 3. In 
spite of the assumed factor of ~ \f?> larger statistical errors the sensitivity of 
SNO to the violation of the equivalence principle is essentially equal to that 
of Super-Kamiokande. It is due to an advantage of the deuterium detector 
in which one can probe neutrino spectrum more accurately than the light- 
water detector by use of the absorption reaction v e + d — > p + p + e~. If we 
ignore the three effects mentioned before the electron energy is equal to the 
neutrino energy minus 1.44MeV. On the other hand the measurement of the 
recoil electron energy implies taking convolution (smearing) of the original 
neutrino spectrum. 

In Fig. 5 and 6 we present the results with the parameter (ii) which cor- 
responds to the large-mixing-angle solution. Fig. 5 is for Super-Kamiokande 
and Fig. 6 is for SNO. Again the results are almost the same. It is evident 
that the sensitivity is at most Af ~ 10~ 14 at some particular values of sin 26q 
and in general | A/| > 10~ 13 . 

What is the reason for such (relative) great insensitivity in this case? 
In what follows we give a very simple-minded argument which apparently 
explains this behavior. For the parameter set (ii) with Af ~ 10~ 14 , the 
adiabatic condition is almost completely satisfied except for the case 6m = 
9g- An angle which diagonalizes the Hamiltonian matrix ( |T3"D is given by 



where Sm,g = sin26 ) A/,G an d Cm,g = cos29 m,g- By using this 6, the two 



sin 26 



AmSm + AgSg 



(25) 
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energy eigenstate V\ and v 2 of H is given by 

y 2 

At the center of the sun, 9 ~ 7r/2 because the matter effect iV e dominates 
over the mass term and the gravitational term (see eq. (|25f)). This means 
that v e produced in the center of the sun is essentially equal to the state v 2 
(see eq. (|26|)). On the other hand, at the earth 6 ~ 6 M because there is 
no matter effect (N e = 0) and A^ <C Am for the parameter set (ii) with 
Af ~ 10 -14 . Hence if the resonance occurs adiabatically the probability of 
observing u e at the earth is given by 

P(u e ) ~ | (u e \u 2 ) \ 2 earth = sin 2 d ear th ~ sin 2 6 M - (27) 

Of course this is the same expression as we obtain in the large mixing angle 
MSW mechanism by which the flux of solar neutrinos is uniformly reduced 
independent of energy. Thus we expect the non-universality of the neutrino 
gravitational couplings is unlikely to produce detectable effect in modulation 
of the neutrino energy spectrum with | A/| smaller than 10 -14 , if adiabaticity 
condition holds. Far better sensitivity with the parameter (i) is then under- 
stood to be due to the non-adiabaticity of the neutrino flavor-gravitational 
transformation. 

If we increase the value of |A/| beyond ~ 10~ 13 the sensitivity to the 
non- universal gravitational coupling of neutrinos depends on the sign of Af. 
As we see in Fig. 5 and 6 it produces an appreciable difference with the 
large-mixing-angle parameter (ii). Let us examine negative Af case first to 
understand the reasons for the difference. In view of the resonance condition 
flIE|) the resonance ceases to occur if Af < and |A/| > 10~ 13 because the 
gravity effect (third term of flTBD ) tends to dominate over the matter effect 
(first term). If the resonance does not occur the spectrum of course deviates 
from the one expected by the MSW mechanism. It explains the region of 
sensitivity which starts to develop at —Af ~ 5 x 10~ 14 in Figs. 5b and 6b. 

In case of positive Af the situation is quite different. If A/>10~ 14 
the resonance can not occur inside the sun because both the first and the 
third terms of ( p^6|) give positive contributions. So the resonance occurs 
between the sun and the earth. For this resonance the adiabaticity condition 
is fulfilled and the argument leading to ( p7| ) holds. The insensitivity exhibited 



cos 6 -suit/ 
sin 6 cos 6 



(26) 
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in Figs. 5a and 6a follows because of the difficulty in perturbing the adiabatic 
large-mixing-angle mechanism by the effect of Af. We should note that our 
computation is not very accurate in regions | A/| > 10 -14 for these particular 
figures because we only integrate the equation (|12|) up to r = 10r Q . 

In the analysis in this paper we have assumed that the MSW mechanism 
is the cause of the solar neutrino deficit. How reliable is this assumption in 
the light of existing data of ongoing four experiments? An extensive analysis 
by Hata and Langacker [FJtJ seems to indicate that the MSW mechanism is 
the most plausible solution to the solar neutrino problem. Of course, we have 
to wait for the next generation experiments, Super-Kamiokande, SNO, and 
BOREXINO, to see if this really is the case. At the present stage of solar 
neutrino experiments it is worthwhile to examine the expected sensitivity 
with other solutions to the solar neutrino problem. We hope that we will be 
able to return to this problem in the near future. 

In summary we have discussed the possibility of testing the equivalence 
principle by using solar neutrinos. We have presented an appropriate for- 
malism for treating neutrino propagation under the influence of the weak 
gravitational fields of the sun. Assuming the MSW mechanism as the origin 
of the observed solar neutrino deficit we analyzed a complex process of neu- 
trino flavor transformation with coexisting effects of the flavor mixing and 
the assumed non-universal gravitational couplings of neutrinos. We have 
obtained, as expected sensitivities, |A/| ~ 10~ 15 — 1CT 16 for a parameter 
corresponding to the nonadiabatic solution and |A/| > 10 -13 for the large- 
mixing angle solution. We emphasize that it opens a remarkable possibility 
that solar neutrino observation by next generation (light- or heavy-) water 
Cherenkov detectors can improve the experimental bound on violation of the 
principle of equivalence by 3-4 orders of magnitude. 
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Figures Captions 



Fig. 1. Modulation of (a) recoil electron energy spectrum and (b) the ratio of 
electron energy spectrum to that of SSM to be seen at the Super-Kamiokande 
detector are indicated. The ratios in (b) are normalized at E e = 9 MeV. 
The flavor mixing parameters are taken as (i) Am 2 = 6 x 10~ 6 eV 2 and 
sin 2 29m = 0.01 as a representative of the nonadiabatic MSW solution to the 
solar neutrino problem. The presented five lines are for Af = 0, ±5.0 x 10~ 16 , 
and ±10~ 15 with corresponding line symbols shown in the figure and 9q is 
taken to be equal to 9 M . The upper solid line in Fig. 1 (a) which is denoted 
as SSM represents the energy spectrum expected by SSM. 

Fig. 2. The same as in Fig. 1 but for SNO. 

Fig. 3. The regions of sensitivity at 90 % confidence level to the violation of 
the equivalence principle are drawn on the parameter plane of Af — sin2# G -. 
The result is for Super-Kamiokande and, 40000 (open circles) and 20000 
(asterisks) solar neutrino events are assumed, roughly corresponding to its 
operation over 4 and 2 years, respectively. The flavor mixing parameter 
for nonadiabatic MSW solution are taken as (i) Am 2 = 6 x 10~ 6 eV 2 and 
sin 2 2# M = 0.01. 

Fig. 4. The same as in Fig. 3 but for SNO. The numbers of events are 
taken as 12000 (open circles) and 6000 (asterisks), roughly corresponding to 
its operation over 4 and 2 years, respectively. 

Fig. 5. The 90 % confidence level regions of sensitivity are drawn for flavor 
mixing parameter (ii) Am 2 = 2x 10~ 5 eV 2 and sin 2 29 M = 0.7 of large- mixing 
angle MSW solution. The result is for Super-Kamiokande and assumes 40000 
(open circles) and 20000 (asterisks) solar neutrino events. 

Fig. 6. The same as in Fig. 5 but for SNO. The numbers of events are taken 
as 12000 (open circles) and 6000 (asterisks). 
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